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ABSTRACT 



The time development of equal-time correlation functions in quantum mechanics and quantum 
field theory is described by an exact evolution equation for generating functionals. This permits 
a comparison between classical and quantum evolution in non-equilibrium systems. 



*^This work was performed in part at CERN, Geneva. 



In a statistical description of nature only expectation values or correlations are observable. 
The dynamical laws describe how a given (complete) set of correlation functions at some initial 
time to has evolved at some later time t. In a very general context this constitutes a differen- 
tial evolution equation for the correlation functions or their generating functionals. Different 
dynamical laws, i.e. the difference between a classical and a quantum description, manifest 
themselves only in different evolution equations. For quantum mechanics of a few degrees of 
freedom we are accustomed to use the Schrodinger equation for pure states and to evaluate 
time-dependent expectation values in an ensemble specified by a density matrix. The classical 
counterpart are Newton's equations with the ensemble described by a probability distribution 
for initial values. For many degrees of freedom a more direct formulation of the dynamics 
in terms of an evolution equation for generating functionals of correlation functions may be 
advantageous. Then one only treats with the relevant observable information and may avoid 
additional complication in intermediate steps of their computation. For classical equations of 
motion such an evolution equation has been established recently 0. It is the purpose of this 
note to develop the counterpart for quantum mechanics and quantum field theory. 

The resulting evolution equation in quantum field theory is a functional differential equation. 
Its practical use depends on the ability to find realistic truncations for the time-dependent 
effective action - the generating functional for the lPI-Green functions. Nevertheless, already 
at the present formal stage a comparison between classical and quantum evolution equations 
sheds some light on current dynamical simulations of non-equilibrium quantum field theory 
problems within a classical approximation for the field equations. It also provides a systematic 
framework to compute quantum corrections to a classical evolution, as relevant, for example, 
for inflationary periods in early cosmology. 

Consider a system with an arbitrary number of degrees of freedom described by conjugate 
operators Qi, P iy i = 1, . . . n, [Q i7 Pj] = iSy, [Q iy Qj] = [Pi, P 3 ] =0, (H = 1). The HamiltonianQ 
is given by 

where P 2 = Pi Pi is a scalar product (summation of repeated indices is always implied) and 
V[Q) = V(Q\ . . . Q n ). In the Heisenberg picture the operators Q and P depend on time 
according to 

Qt = -Pi, Pi= -Vi[Q] (2) 
m 

(We use a notation Vi = J^, = qq.qq , etc.). Let us now define the generating functional 

Z[j,h,t] = Tr(e jQ(i)+/lPW p) (3) 

where p is the density matrix which is time-independent in the Heisenberg picture. It is a 
function (or functional for n — > oo) of the time-independent sources ji and hi. The coefficients 

1 Systems with different masses mi can be brought to this form by simultaneous rescaling of Qi and Pi - in 
fact, we could also use a scaling with m = 1. Our formalism can be generalized for arbitrary H[P, Q] but we 
restrict the discussion here to the most common form quadratic in the momenta. 
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of its Taylor expansion 
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(4) 



are time-dependent symmetrized expectation values of correlation functions. For a one-component 
example (n = 1) 

z (M) (t) = <(Q*i*) s (t)> = Tr {(Q k (t)P\t)) s p} (5) 

the symmetrized ordering ( ) s stands for an equally weighted sum over all possibilities 
to form different chains of k operators Q and I operators P, i.e., 



(Q 2 P 2 ) S = \{Q 2 P 2 + QP 2 Q + P 2 Q 2 + PQ 2 P + QPQP + PQPg) 
o 



(6) 



The generalization to arbitrary n is obvious. Knowledge of Z[j,h,t] contains the complete 
information about the system. Macroscopic quantities and thermodynamic functions can be 
expressed in terms of correlation functions 0| . 

The time- dependence of Z obeys the evolution equation (d t Z is the time derivative at fixed 
j and h) 

d t Z = tTr([H,e^ +hP } p) (7) 

Our aim is to express the right-hand side in terms of Z and its derivatives with respect to j 
and h. In order to deal with the problem of ordering the non-commuting operators we use 



i[H,e^ +hP ] = - ]t -?-e^ hP + iV 
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m dh 

and evaluate the commutator with the Campbell-Barker-Hausdorff formula 

e /3Q ; e jQ+hP _ ( e iph _ e -|/9fc^ e {P+j)Q+jP (g) 

For each of the two contributions on the right-hand side, one can express ^ by an appropriate 



combination of 4^ and h. This leads to 

9] 



V[Q],e 



jQ+hP 



i V 



d i , 
Sl + 2 h 



V 



d i , 
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(10) 



and (H), ( |10|) can now be inserted in the trace (0). Our final result is a linear partial differential 
equation for the evolution of Z with a simple structure 



d t Z = c {z) z 

C (z) 
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For a comparison with the time evolution of Z in classical statistics JT|] we expand V in 
powers of h 
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(12) 



The first term in this expansion reproduces the classical evolution equation. This leads to the 
important observation that for linear equations of motion (H quadratic in Q and P) there 
is no difference in the time evolution of correlation functions between quantum statistics and 
classical statistics! The higher order terms appearing for non-linear equations of motion can 
be viewed as quantum corrections to the classical evolution. Restoring h, they involve powers 
h 2 , ft 4 , etc. For V containing only up to quartic terms the quantum correction ACqIj reads 
explicitly (V ijk = V ijk [0]) 
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(13) 



It may be instructive to consider two examples. We first take a single anharmonic oscillator 
with 

i i 
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V(Q) = - muj 2 (Q 2 + -V^mJ^Q 3 + 2mu5Q A 



In terms of a complex source 



J 



1 



--(j + imujh) 



(14) 



(15) 



ylmuj 

which is conjugate to the creation operator of the harmonic oscillator the evolution equation 
reads 



( d d\ J d d 



d.J* d.J 



dJ* dJ . 



(16) 



In our picture the state of the system at a given time t is described by Z[J, t]. For the harmonic 
oscillator (7 = 5 = 0) the stationary states or fixed points of Z (d t Z = 0) are exactly those 
for which all terms in Z involve an equal number of powers of J and J*. They correspond 
to incoherent mixtures of eigenstates of the number operator a^a or energy eigenstates. It 
is intriguing that due to the equivalence of the quantum and classical evolution for linear 
equations of motion the quantum mechanical energy eigenstates can also be viewed as classically 
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stationary probability distributions for commuting coordinate and momentum variables. For 
7,5^0 these particular states are not stationary any more. We know, nevertheless, that 
Eq. (fl6|) must admit an infinite number of fixed point solutions which correspond exactly to 
incoherent mixtures of energy eigenstates of the anharmonic oscillator. This follows ^from the 
general observation that Eq. ([/]) is equivalent to 

d t Z = -*Tr {e jQ+hP [H,p}) (17) 

and H(t) = H(to). Every p which commutes with H defines a stationary state. An important 
qualitative difference between the quantum and classical evolution equation can be easily seen 
if we neglect 5. The quantum mechanical evolution equation for InZ contains a "correction" 
term jo(J* — J) 3 which acts as an additional constant "force" on the connected correlation 
function for P 3 , i.e., Ad t < P 3 > c = ^u! r y(2mco) 3 ^ 2 . A cubic potential is not bounded from 
below. The local minimum at the origin is separated from the unstable part by a barrier. 
A classically stable initial ensemble with all energies below the barrier becomes unstable in 
quantum mechanics due to tunneling, as reflected by the additional "force". The essential 
features of this effect are not changed if we restore stability by a small positive non-zero 5. 
Once reexpressed in terms of real sources j and h, Eq. (|1^) is a real partial linear differential 
equation for a function of three variables Z(j,h,t). A linear combination of two solutions is 
again a solution. A numerical solution of this equation contains at once the information about 
the time evolution of all expectation values of arbitrary powers of Q and P\ In particular, it 
seems interesting to investigate the outcome for classically chaotic systems. 

Our second example is a linear chain of oscillators (with mass m scaled to one and d = 1) 

H = H{l p i+ \fQ\ + \™ d/2 Ql + \^~ d Qi + ^ - Qi) 2 } (is) 

Here a is the distance between two oscillators on the chain and we take periodic boundary 
conditions with a fixed length Q = Yh=i The evolution equation reads 




Again, the last two terms are the quantum corrections. 

The transition to a field theory can be made by taking the limit a — > while keeping a 
non-zero "volume" Vt. With the replacements a~ d ^ 2 Qi —> <p(x), aT d l 2 Pi = tt(x), a d J2i ~^ 
J d d x, a~( d ' 2+1 ) (Qi+i — Qi) — > d(p/dxi - we generalize here to an arbitrary dimension d - the 
Hamiltonian fll8|) becomes 

H = j d d x^-^{x)+ l -d^{x)d^{x) + U{^{x))^ 
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urn*)) = ^ 2 w + ^ 3 w + ^ 4 (x) 



(20) 



This is the Hamiltonian of a relativistic scalar field theory with field equations 7r(x) = tp(x) 
7f = <p = A(p — dU /dip. Since also the commutation relation (a~ d 8ij — > 8(x — x')) 



[<p(x),tt(x')] = i6(x — x' 



(21) 



is covariant under Lorentz transformations, our system describes a Lorentz-invariant quantum 
field theory. Here Z becomes a functional of the sources j v (x), j n (x) which are related to the 



sources in the discrete version by a 
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Tr{exp(J d d x[j lfi (x)(p(x) + j n (x)7r(x 



dhi Sj n (x) ' 

time evolution equation for a scalar quantum field theory 



dji Sj v (x) 

We finally arrive at the 



d t Z[j v ,j n ,t] 



d d x\j v (x) 
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Z[ Jv , Jn ,t} (22) 



Its generalization to several scalar fields <p a (x) is straightforward - with tt 2 — > ifa^a, di(pdi(p — ► 
di(p a di(p a in H (|TJ) and commutation relation [<p a (x), ftb(x')] = i8(x — x')8 a b we only have to 
replace j v -±- -> j^^, J^A^- -> J^A^- in Eq. (^). All symmetries oiU(<p) are preserved 
by the evolution equation in the sense that an initially symmetric state remains so at later 
time. This does not preclude spontaneous symmetry breaking in the course of the evolution - 
it may be detected by adding a small symmetry breaking linear term in U or by starting with 
a slightly asymmetric initial state. 



Equation (2"2"|) is a functional differential equation and its approximate solution has to pro- 
ceed by some truncation. For this purpose it seems an advantage to switch to the generating 
functional for the lPI-Green functions, F[<p, 7T,t] = — In Z[j,t] + J d d x(j (p (x)ip(x) + j n (x)ii(x)), 
where (p(x) = Jp^Jy, ir(x) = J^rfj ■ The derivation of the evolution equation for the effective 
action T proceeds as in Ref. [|TJ and we only give here the result for the Hamiltonian ( p0|) with 
d t T a time derivative at fixed tp and 7r 



d t T 
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v(<p (x) + G w (x, x)) + - lip (x) + 3ip(x)G ipv> (x, x) 
dx\dx 2 dx^G tsrix \ X - X\)G ^2 (•£; 2^2) 
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Here G 7 y(x, x') = (<p^(x)(py(x'))j — <^ 7 (x)<^ 7 /(x') is the propagator in presence of sources where 
<y3 7 , 7 = tp, 7r is a shorthand for (<£>, it). The propagator can in turn be expressed by the inverse of 
the matrix0 of second functional derivatives of T, i.e., G rr r(x, x') = (r( 2 ))~y(x, x'). We observe 

that Cj) plays the role of the classical Liouville operator, with (p n replaced by (<fi n ). The 
difference {(p n ) — ip n is accounted for by the terms involving G. They induce a dependence of £ 
on T and turn the evolution equation non-linear. The quantum corrections are all proportional 

to (g ) 3 - 

The evolution equation (|23]) has a fixed point r*(/3) corresponding to thermal equilibrium 
where p = Z^e-? 11 , Z = Tre- pH in Eq. (§. This can be computed by functional integral 
methods ||. It is far from obvious, however, if and in what sense the solutions of Eq. (|23| ) 
with non-equilibrium initial conditions approach this fixed point. A uniform approach is not 
possible due to the existence of infinitely many other fixed points which correspond to incoherent 
mixtures of eigenstates of H . At best, the equilibrium fixed point can be approached if we 
restrict the discussion to correlation functions of suitably averaged fields (coarse graining) or 
to other subsystems. Thermalization can also be achieved by a coupling to an environment, 
thus introducing a stochastic element in the equations of motion. A truncation of T[ip, ir,t] 
may destroy the existence of the infinitely many fixed points. Information about higher 1PI 
correlations or their precise momentum dependence is omitted in this way. It is conceivable 
that truncated equations have a more uniform approach to the fixed point than the exact 
ones. Good truncations should at least retain those terms that play an important role in the 
computation of r*(/3) by the solution of renormalization group equations in dependence on a 
coarse graining scale Q. A better understanding of the impact of truncations will be crucial 
for the practical use of the present formalism. 
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